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Abstract 

This paper continues the preceding paper on the problem of quantum dynamics on the lattice. 
Firstly we consider the multiple reflections of the wave function (Loschmidt echo) . The phenomenon 
of wave function concentration on the impurity site after reflections is found. The solution rep- 
resenting the total amplitude a{t) is obtained as the series in terms of partial amplitudes ak{t). 
The contribution of kth. partial amplitude becomes dominant only after kih. reflection from the 
lattice end. An excellent agreement between analytical and accurate numerical results is obtained. 
Next problem, - wave packet trapping by defects, is solved by numerical simulation. Analytical 
expressions are derived in few cases allowing to estimate the quantum efficiency of charge trans- 
fer. Obtained results can qualitatively explain recent experiments on the highly efficient charge 
transport in olygonucleotides and polypeptides. 
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I. SETTING UP A PROBLEM 



For the completeness we repeat few details of the setting up the problem from the pre- 
ceding paper, but in a somewhat different notations. 

We consider quantum system consisting oi {N + 1) sites. Most left site is an impurity 
site; other N sites are reservoir. This system is defined by the (A^ + 1) x (iV + 1) matrix 
hamiltonian 

where E is the on-site energy of the impurity site; B is the N x N tridiagonal matrix of 
reservoir with matrix elements B^j = -|- 5i,i-i; ^^ is the interaction of the impurity 
site with the reservoir and is chosen in the form vi = CSi^i. The wave function is ^(t) = 
|a(t), l)(t)| = a{t),bi{t),b2{t),...,bN{t). Initial condition: a{t = 0) = 1, '^{t = 0) = 0. 
The aim of this paper is to find the amplitude a{t) on the impurity site as the result of 
multiple pave packet reflections (Loschmidt echo). 
The dimensionless {h — 1) Schrodinger equation is 

i^^^Ea{t) + {'^ -^y, a(i = 0) = l 

4^ (2) 



dt 



(B-~S{t)^ +^ a{ty, b(i = 0) = 0. 



If vector b is expanded in terms of the cigcnfunctions bk{i) of matrix B, then the following 
integro-differential equation can be obtained: 

t 

a = -iEa - J B^{t - r) a{T) dr, a{t = 0) = 1, (3) 



with the kernel 

^ — > 2 

5^(t) = J]expH£,t)(^-bk) (4) 

k=l 

and £fe is kth eigenvalue. 
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The Schrodinger equations 



for this model are: 
ia = Ea + Cbi 
ihi =1)2 + Ca 
ib2 =61 + 63 



(5) 



ibN = bN-i 

with initial conditions a(0) = 1, 6j=i 2,...,Af(0) = 0. Eigenf unctions e{k) and eigenvalues 6j(/c) 



of the matrix B are well known: 

irk 



e{k) = 2 cos 



+ 1 



6fe(i) 



+ 1 



sm 



k 



N + 1 



k = l,2,...,N. 



(6) 



The kernel B (t) of Eq. ([3]), in accordance with ([6]), is defined by the following sum: 



2C' 

WT 



N 



k=l 



sm 



7ck 



N+1 



exp 



— 2i cos 



Tck 



N+1 



(7) 



II. EXPANSION OF a{t) IN TERMS OF PARTIAL AMPLITUDES 

In the preceding paper we have shown that the well formed impulse with sharp forward 
front is generated at large times. The velocity of the impulse front f = 2 is the maximal 
group velocity. If a lattice consists of sites, then the time of the first returning to the 
impurity site is t ^ N. And the cycle of k returnings takes t ~ kN. 

In this paper the detailed analysis of amplitude a{t) in the time range t ^ N was done. 
Now our goal is to find amplitude a{t) on the impurity site at all times. We solve equation 
(j3]) with the kernel given by ([7]). 

Bearing this in mind, the solution of (|3]) is searched as the expansion in terms of partial 
amplitudes ak{t). Each amplitude afc(t) is negligible in the time range t < kN and its 
contribution to the total amplitude a{t) becomes essential starting only from times t ~ kN. 

In order to make sure in this possibility, the kernel B^[t) of Eq. ([7]) should be transformed 
according to the Poisson formula: 

m=+oo 



B^{t)=C^ Bm{t), where 

m=—oo 

TV 

Emit) = ^y" dxsin^(x)exp{-2i[m(A^ + 1) +tcos(x)]}. 



(8) 



Note that Bm{t) is the real function and there exists the exphcit expression for -Bo(^) 
through the Bessel functions: 



Bo(t) = Jo(2t) + J2(2t). (9) 
The original equation (|3]) with regard to transformation ([8]) is rewritten in the form: 



t 

oo 



a{t) = -iEa{t) -C^ / ~ 

m=— oo Q 

Note an important and useful property of functions Bm'- every function B^it) is small in 
the time range < t < mN, such that 

, , sin(2t) , , , ^ 

With a reasonable degree of accuracy one can say that the function Bm{t) ~ when 
t < mN. Consequently, when t < mN, then only terms Bk with indices k < m are essential 
in the governing equation ( ITOl) . 

Now we represent amphtude a(t) as the sum of partial amplitudes 

a{t) = ao{t) + ai{t) + a2it) + . . . . (12) 

Recall that our assumption concerning amplitudes ak{t) lies in the fact that amplitude ak{t) 
is small at t < kN. Their contribution to the total amplitude a{t) becomes essential starting 
from time t ~ kN. To some extent afc(t) behaves like Bk{t). Deriving an equation for ak{t), 
we assume that its right-hand side contains only terms Bm{t) with m < k. 

Prior to write down an equation for ak{t), we make one comment. Terms Bm(t) are small 
(even for not too long lattices). This smallness is ensured by the fact that B^rnif) = B„i{—t). 
And B.fn{t) ~ (mN)^^ in the considered time range {t > 0). Nevertheless these terms are 
taken into account as their accounting allows to derive explicit analytical formulaes. 

Lets introduce following functions to account the contributions from Bm{t) (m < 0): 

B^{t) = Brn{t) + 5_^(t), m > (13) 

and for uniformity of notations we set -Bo(^) = B^lt). 

Now in accordance with our assumptions, convert one equation ( ITOj) in a system of cycling 
equations for the partial amplitudes ak{t). Equation for ao(t) contains in the right-hand side 



only function -Bo(^) (see ([9]) 



do = -iEao - j Bo{t - t') ao(t') dt', ao(0) = 1. (14) 



and in accordance with the above-mentioned arguments, an equation for ak{t) {k > 0) can 
be rewritten: 

ak = -itjau - / B^{t - t') ak-m{t') dt', 0^(0) = (A; > 0). (15) 



lEau -C^ ^-(^ - (0 dt', afc(O) = (A; > 0). 

Q m=0 



It is notable that system fllSp is accurate. 

The Laplace transformation of f[T^ gives the following expression for ao(p)- 

1 



ao{p) 



p + tE + C^Bo{p) 
Here Bq{p) is the Laplace transform of -Bo(^) (see Q): 

Boip) = I (v7+4 



(16) 



p 



(17) 



Recall that expression (fT6l) is nothing else then amplitude a{t) in the infinite lattice when 
there is no reflections. 

The Laplace transformation of system (fT5|) gives the algebraic system of recurrence rela- 
tionships for ak{p) (A; > 0): 



fc-i 



(p) = -C^ ao(p) ^ Bk^rn{p) ttmip)- 



(18) 



m=0 



The Laplace transform Bm{p) of the function B^it) (see (|8])) is given by the integral: 

t 



7r 



Bniip) = - / sin^(?/) exp[-2im(A^ + 1)?/] 



1 



+ 



1 



p + 2i cos{y) p — 2i cos{y) 



dy. (19) 



An explicit form for the integral f lTI?]) can be obtained. With this aim in view, the inte- 
gration contour should be deformed as shown in Fig. [TJ As the integrand in fll9p has period 
71, integrals Ji and I2 are cancelled as they have opposite integration paths. Consequently 
Bm{p) is defined only by the pole contribution at the point y^ (cos(?/o) = w/'^)- Then we 
have: 

_ r _ -1 2m(Af+l) 

Bra{p) = V4+^ iBoip) . (20) 
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FIG. 1: The integration along the path [0, tt] in expression (jl9p is changed to the residue in the 
pole po and two integrals Ii and l2- 

Due to the fact that B^ip) forms the geometrical progression by m, it is possibile to make an 
explicit summation in the recurrence formulaes ( !T8|) . It is the reason why the total amplitude 
a{t) can be represented in the form of the closed expression (this is done in Appendix). 
The Laplace transforms of the partial amplitudes ak{p) {k > 0) are expressed as 

2fc(iV+l) 



= -ao{p)d{l - d) 



k-l 



iBr 



where d = uqC'^ \/p'^ + A. 



(21) 



To get explicit expressions for the partial amplitudes, the inverse Laplace transform 
should be made. Transforming the Laplace integral to the integration path around the cut 
[—2i,2i], one can get an expression for ao(t): 

2 



ao(t) 



IT 



exp{ixt) Im 



X 1 



E 



dx. 



(22) 



Further on we get the following expression for the partial amplitudes akit): 
2 

an I i — 



(J2 

ak{t) = / "\/4 — exp{ixt)Re 



TT 



1 - C^aaVA 



X' 



k-l 



iBc 



2k(N+l) 



dx. (23) 



Here the Laplace transform for ao(p) (see f|T6|) and fl20|) ) are performed sX p = ix. An 
analysis of expression for afc(t) shows that, as supposed, afc(t) ^ 1 when t < {k + 1)N. 

Therefor, if the partial sum ao + ai + ... ak is taken for the representation of amplitude 
a{t), then the error of such approximation is of the same order as the smallness of Bk{t) 
(see ([II])), ~ [k + 1){N + 
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FIG. 2: The comparison of the hmited sum of partial amphtudes + ai + a2 (solid line) with 
the nmnerical integration (dashed line). Dotted lines ~ partial amplitudes oq, ai, 0,2 "starting" at 
times t = 0, 2, 4, correspondingly. Small divergence is observe only at t > 8 where the unaccounted 
partial amplitude 04 (dash-dot line) starts to make the contribution. Parameters: N = 2, E = 
1, C = 0.5. Mean-square error (MSE) < 10^^. 

Consider as an example very short lattice (A^ = 2), and as an approximation - sum of 
only three partial amplitudes + cii + 02- In Fig- EJthis partial sum is compared with the 
result of numerical integration of the Schrodinger equation ([5]). For t < 8 the expected error 
is ~ 5 ■ 10~^. Thus the representation of a{t) by the sum of partial amplitudes is a very good 
approximation even for short lattices (an accuracy increases if lattice is longer). 

III. RECURSION. MULTIPLE RETURNING TO THE INITIAL STATE. 

If the lattice is long enough then partial amplitudes, following each other, have enough 
time to damp on the corresponding time range [t-^t + N]. In this case the partial amplitudes 
do not interfere and reproduce the total amplitude with very high accuracy (see Fig. |3]). The 
maxima of returning amplitudes slowly decrease. 

Partial amplitudes interfere in the short lattices and maxima of returning amplitudes are 
irregular. The dependence of the total amplitude a{t) vs. time for the lattice with = 10 
is shown in Fig. HI 
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FIG. 3: Sum of partial amplitudes oq, oi, 02, 03 practically coincide with the total amplitude a{t). 
Numerical results are not shown as they excellently coincide with analytical result (MSE < 10"^). 
Parameters: TV = 100, E = 0,C'^ = 0.25. 




FIG. 4: Solid line - sum of partial amplitudes |ao + Qi + • • • , +o.8|- It practically coincides with 
amplitude a{t). Dots - partial amplitudes ao, ai,...,08. Maximal value of returned amplitude 
is 0.95 (at t = 42). Main contributions to maximum give partial amplitudes 02, 03, 04, as, ae. 
Empty circles - numerical result. MSE < 10~^. Parameters: N = 10, E = 0, C'^ = 0.4. 
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Numerical analysis performed at different values of parameters C, E N shows, that the 
maximal value of returned amplitude Oret ~ 0.972 at t = 42 for N = 10, E = 0, C'^ = 0.4, 

Incident and reflected impulses interfere on short lattices and the degree of returning is 
difficult to analyze at arbitrary parameter values C, E, N. But the first returning (maximal 
value of the partial amplitude Oi) can be treated analytically if the lattice is long enough 
when amplitude oq becomes negligible. 

The expression fl23|) for ai{t) using the trigonometric substitution of variables can be 
written: 

ai{t) = / sin^ X exp 2zt cos X Re ^ / \. dx. 24 

vr J [ [2 cos(x) + i? — C"' exp(2x)J'' J 



Fig. [5] shows the maximal values of amplitude calculated according to (!24l) at E = and 
different values of parameters and C. One can see that if ^ 0.2 then the returned 
amplitude practically does not depend on the lattice length (10 < < 100). The dissim- 
ilarity of the partial amplitude from the total amplitude is negligible on this time range 
(A^ < t < 2N). Divergence becomes essential (~15%) for the shortest of considered lattices 
(A^ = 10) and smallest value of parameter C (C^ = 0.1). Amplitude Oq has no enough time 
to fully decay at these parameters values. 

IV. WAVE PACKET TRAPPING BY AN ACCEPTOR 

The phenomenon of multiple reflections of the wave packet (Loschmidt echo) is unlikely 
to observe experimentally. The reason is that the wave function does not interact with 
environment. Below we consider the problem which mimics the experiments on the charge 
transfer (CN) in DNA where the wave function is irreversibly trapped by an acceptor. And 
the fraction of the wave function trapped by an acceptor is of primary interest of this section. 
This quantity can be compared with the quantum efficiency of CT. 

Consider the lattice with the attached site (acceptor, see Fig. [H]). The number of this 
site is A'a. The acceptor has on-site energy i?a and the hopping integral Ca. Amplitude of 
the wave function on the acceptor is labelled by 6a- Initially we limit ourself by the weak 
bounding energy between the lattice and acceptor, i.e. Ca ^ 1. This approximation allows 
to make necessary analytical estimations. 
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FIG. 5: Maximal value of amplitude of the first returning at E = and different values of C and 
N. 




a 



FIG. 6: Schematic representation of the lattice with acceptor. Acceptor is attached to the lattice 
site with number and has the on-site energy Eg, and hopping integral (the interaction energy 
with the lattice) Ca. 

The system of equations ([5]) changes in an obvious way: an equation for the amplitude 
of the wave function on the acceptor is added: 



Kit = 0) = 0. 



(25) 



The equation for the site N^^ is also modified: 



(26) 



Other equations stay unchanged. 
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FIG. 7: The dependence of the wave function amphtude on the acceptor |6a(t)| for two positions 
of the acceptor on the lattice: A'^a = 5 (sohd hne) and A'^a = 15 (dots). Dashed hne - expression 
(pl). Parameters: N = 100, C = 0.4, E = 0.5, Ca = 0.02, = 0.3. 

Staying in the frameworks of the initially formulated problem, consider now results ob- 
tained in numerical simulation in the case, when the acceptor is located close to the impurity 
site (donor). It turns our that the fraction of the wave function on the acceptor, being cap- 
tured, stays on the acceptor for a long time. Fig. [7] shows this phenomenon for the lattice 
consisting of = 100 sites. (The time range is such, that the reflected impulse has no time 
to return back after reflection). 

It is possible to estimate the dependence of the wave function amplitude on the acceptor 
vs. time in the approximation of the weak coupling (Ca ^ 1). Note that in this approxima- 
tion the acceptor affects the lattice very weakly. Therefor the lattice is not disturbed and it 
is described by Eq. ([5]). Amplitude of the wave function on the acceptor will be calculated 
according to the obvious expression resulting from fl25l) : 





Fig. [8] shows the comparison of two solutions of the Schrodinger equation: accurate 
(expression ([5]) ) and approximate (equation (125|) ). One can see that these solutions differ 
very little and an approximation by the unperturbed lattice is very good. 




(27) 
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FIG. 8: A comparison of the accurate and approximate solutions for the acceptor attached to 
the tenth lattice site (A'^a = 10). Solid line - accurate solution (an influence of the acceptor on the 
lattice is accounted). Dashed line - unperturbed lattice. Parameters: N = 100, N^, = 10, C = 
0.4, E = 0.5, Ca = 0.02, = 0.3. 

Analytically will be considered the case when time is large enough such that amplitude 
a{t) on the impurity site decreased practically to zero. It means that the lattice is long and 
time is large, and the impulse and its tail went away from the acceptor. Then the lattice 
can be considered as having infinite length. 

As is seen from (127|) . it is necessary to evaluate integral jQexp{iE^T)bN^dT (the phase 
multiplier exp{—iEi^t) is unessential for the modulus of the wave function). We consider the 
amplitude on acceptor in the limit t — )■ oo. 

To evaluate the integrals, system ([5]) should be multiplied by exp{iEt) and integrated in 
the limits from to oo. Lets introduce the notations: 



oo 



oo 




(28) 
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Then for Ik we get the recurrence relations: 

Ch=-t + iE~ E,)Io 
h = Es^Ii — CIq 

h = EJ2 - Ji (29) 
h = — I2 



Amphtude a{t) on the impurity site in the considered approximation is ao{t) and integral 
Jo is the Laplace transform ao{p) {p = iEg). As the result we get: 

/o = [i[E - E^) + exp(z0)] "\ <p = arcsin(Ea/2). (30) 

System of equations (!29|) has the following solution: 

4 = C/oHexp(z0)]\ (31) 

Thus the limiting values of the acceptor amplitudes (with the accuracy of oscillating multi- 
plier exp(— ii^t)) on different sites are— zCa/fc (see fl26l) ) and differ only by phase multiplier. 
In this case the amplitude on acceptor is 

b^it) ^ -iC^Cexp{-iEj)Io exp[-^exp(^0)]^ (32) 

If the phase multipliers, unessential for the amplitude of the wave function, are eliminated, 
then the amplitude is 

''■W" .(g-E.)%expW) - * = (33) 

Fig. [7] demonstrates that the limiting value of amplitude coincides with the numerical 
calculation. 

It follows from ( l33l) that at fixed values of the hopping integrals C and Ca, maximal value 
of 6a is achieved in "resonance" values of E and i?a, when E = E^. In this resonance case 

&a(t) ^ CJC. 

If times are such that the impulse reflects, returns and passes by the acceptor, then 
the amplitude variations are irregular and depend on the acceptor location on the chain 
(see FigJH]). 
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FIG. 9: Amplitude of the wave function on acceptor when A'a = 5 (sohd hue) and = 10 
(dashed Une). Time is such that the double reflection occurs. Parameters: = 100, C = 0.4, E = 
0.5, Ca = 0.02, = 0.3. 

There was considered above the cases, when the bounding of the acceptor with the lattice 
is weak. But practically, for the efficient charge transfer, it is necessary to obtain the 
conditions when the degree of the CT is higher, i.e. parameter Ca should be larger. The 
value Eg, also plays some role. 

Below we consider few examples when parameter Ca has comparatively large value. 
Fig. [10] shows the dependence of the wave function amplitude on the acceptor, when the 
hopping integral Ca = 0.1. The amplitude becomes well larger and reaches value |6a| ^ 0.2. 

In the case of the total resonance (when E = E^. and C = Ca), amplitude |6a| becomes 
even more and this case is shown in Fig. [HI 

From two latter figures it follows that the acceptor population can change significantly 
depending on the parameters. When the population is small (Fig. [TOj) then time of life is 
comparatively large. And on contrary, time of life is small when population is large (Fig. [TTI) . 
This peculiarity has natural explanation: the larger is the acceptor bounding with the lattice 
the shorter is time of life. 

And finally we consider the case when an acceptor is located on the right lattice end (A^th 
site is an acceptor) (see Fig. [T2|) . It is seen that the amplitude on acceptor is rather large 
(|6a| ~ 0.35), and, what is very important, does not depend on the lattice length. Moreover, 
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FIG. 10: Dependence of the wave function amplitude on the acceptor vs. time for few locations of 
the acceptor on the lattice. A'a = 2 (dashed line), A^a = 5 (solid line) A'^a = 15 (dots). Parameters: 
N = 100, C = 0.5, E = 0.5, Ca = 0.1, = 0.5. (According to §^ the value = 0.2). 









^^^^ ***••••** 
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FIG. 11: Dependence of the wave function amplitudes vs. time on the acceptor |6a(i)| (solid line), 
on the impurity site |a(t)| on the left lattice end (dots) and on the acceptor site |6io(i)| (dashed 
line). Parameters: = 100, C = 0.1, E = 0.5, Ca = 0.1, E^ = 0.5. 

time dependencies for the amplitude decay are practically identical. This result is in good 
agreement with experiments on the charge transfer in synthetic DNA and polypeptides where 
the CT probability does not depend on the distance. 
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FIG. 12: Dependence of the wave function amplitude on the acceptor (located on the right lattice 
end) vs. time for two lattice lengths: N = 100 (solid line) and = 50 (dashed line). The 
time point of reference is shifted back by the value N/2 for the data comparison. Parameters: 
C = 0.5, E = 0.5, Ca = 0.1, = 0.5. 

An estimation of typical time scale is necessary. It has to be done to understand how 
long the wave function stays in the bounded state, and is this time enough for photophysical 
or electrochemical response for the charge registration. The typical dynamical time (period 
of one vibration) is [t]d ~ 1-7 ■ 10"^'^ s jl|, Q]. The typical electronic time (time unit in this 
work) is approximately two orders of magnitude shorter [t]e ~ 2.2 ■ 10^^^ s j^. As is shown 
above, the wave function can stay on the acceptor during dozens time units, what is ~ps. In 
many cases this time is enough for effective charge trapping by an acceptor with following 
registration. 



V. CONCLUSIONS 



In two papers we thoroughly analyzed the quantum dynamics of the excitation (electronic 
wave function) propagation (first part), reflection and trapping (second part). The system 
consists of the homogeneous one-dimensional lattice with the impurity site, and an excitation 
initially is totally localized on the impurity site. 

A rather unexpected results consists in the fact that initially localized wave function starts 
to move spontaneously forming well defined wave packet. After first refiection wave packet 
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is again concentrated on the impurity site with the amphtude > 90% of initial amphtude. 
This process repeats many times. 

To describe multiple reflections of the wave packet an useful approach consisting in the 
representation of the full wave function on the impurity site a{t) through the partial ampli- 
tudes afc(t). 

The temporal evolution of the wave function is described with very high accuracy up to 
dozens reflection. The interference of falling and reflected impulses occurs after these large 
times, which is difficult to take into account analytically. The behavior of the quantum 
dynamical system is regular in this time range. The behavior on large times needs further 
detailed consideration. 

Results on the wave function trapping by an acceptor can explain recent results on the 
efficient ballistic charge transport in synthetic DNA and polypeptides. 

Appendix A 

Original equation for the amplitude a(t) has the form (see ([3])): 




(Al) 







The solution of this equation for the Laplace transform a{p) is: 



p + iE + B^ (p) ■ 

The Poisson representation for the kernel B^{p) has form (see ([8]), (|T7j) . (120|) ): 



(A2) 



CXI 



B^ip) 



C^J^BUp) 



m=0 



Boip) 



- (^y/p^+A-p^ 



(A3) 



Bm{p) 



m > 



where the notation 




(A4) 



is introduced. As values Bm form the geometrical progression, we get: 




(A5) 
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And the final expression for tlie total amplitude is: 



a{p) = ao 



1 - b{p) 



(A6) 



l + b{p) (c^ao{p)^/fTi-l 



) 



An expansion into series in terms by b{p) gives Laplace transforms of partial amplitudes. 
By this means it was not necessary to construct the system of the recurrence relations for the 
partial amplitudes, but simply use the expansion a{p) (see flA6|) ) into series in terms by b{p). 
But this approach is valid only for the considered model where the necessary summation 
is possible and the compact expression for the kernel can be written. The method 
suggested in the paper of the expansion into series by partial amplitudes can be applied in 
other problems. 

The back Laplace transformation gives the desired expression for a{t): 



Numerically it was verified that expression (1A6P gives accurate results. 

It worth noting one intriguing property. If the integral (]A7|) is closed around the cut 
[—2i, 2i], then the result (difference of integrals taken along banks of the cut) is zero. It was 
found numerically. And it follows that the function a{p) has poles in the complex plane. And 
amplitude a(t) can be obtained as the sum of residues in these poles. Additional analysis is 
necessary to throw light on this fact. 
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